Comment on "Casimir energies with finite-width mirrors" 
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We comment on a recent publication 1] on Casimir energies for material slabs ('finite width 
mirrors') and report a discrepancy between results obtained there for a single mirror and some 
previous calculations. We provide a simple consistency check which proves that the method used in 
[1] is not reliable when applied to approximations of piecewise constant profile of the mirror. 

We also present an alternative method for calculation of the Casimir energy in such systems based 
on 0]. Our results coincide both with perturbation theory and with some older Q and more recent 
calculations, but differ from those of 
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In a recent publication on Casimir energies and inter- 
action of material slabs ('finite width mirrors') [ij there 
was presented a thorough and rather general treatment 
of the problem at hand. To model the presence of matter 
in the system of massless scalar field the authors of ill] in- 
troduce into the action an additional term concentrated 
in a given domain of the space-time. In the framework of 
renormalizable quantum field theory this method was ini- 
tially proposed by Symanzik ^ though his paper is not 
cited in Such models with position dependent mass 
terms (also called defects) were studied extensively in re- 
cent years, in particular in the case of delta-potentials 
which effectively describe thin films, sec for instance 6], 
. Surprisingly the problem of interaction of finite width 
mirrors hasn't been investigated in full, despite several 
attempts e.g. Q, Hi, [1]. 

The Casimir energy for massive scalar field interacting 
with a homogenous plane slab of finite width has been 
calculated quite a lon^ time ago by Bordag 3 , and also 
recently rederived in |4| by Vassilevich and Konoplv aJlSj . 
A similar problem is considered in Section III.B of [l|. A 
smooth approximation of piecewise constant potential is 
investigated there and Casimir energy for the case when 
this limit is approached is obtained. However the out- 
come presented there (equation (68)) disagree with (well 
defined) massless limit of the above mentioned results, 
and in [l| there is no any discussion of this discrepancy. 

In we elaborated a calculation approach for singu- 
lar potentials describing the interaction of thin material 
films with fields of quantum electrodynamics. In order to 
generalize it to the case of 3-dimensional defects we stud- 
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ied recently the model of massive scalar field interacting 
with a single plane slab . We developed a simple calcu- 
lation method that differs both from [sl, and also from 
one employed in We calculated the propagator and 
the Casimir energy Q . In appropriate regularization the 
later result coincides with ones calculated in Si (but 
differs from that of pj]). It is also in agrement with usual 
perturbation theory both for the massive case and in the 
massless limit. 

One could argue that the discrepancy between results 
of i i i and [l| originates from the smooth approxi- 
mation used in [I] to describe singular (non-continuous) 
potential. Similar problems have been discussed in 
the framework of Dirac equation perturbed by a delta- 
function potential (e.g. see [i3l). However, in our case 
it is not the reason of the above mentioned discrepancy. 
As we show below, the calculation method proposed in 
[l| is not reliable when applied to the case of (approx- 
imated) piecewise constant potential, and the final ex- 
pression (68) [ij is to be reconsidered. 

In Section I of this Comment we present a consistency 
check of the calculation method developed in [l| , while in 
Section II we sketch our own approach to similar prob- 
lem. 



I. CONSISTENCY CHECK 

The calculation approach of [1] is based on expressing 
the Casimir energy as a trace of a logarithm of an in- 
tegral operator, and this step is widely used in Casimir 
calculations. The trace of the operator can be calculated 
directly as a sum of its eigenvalues. For their determi- 
nation authors of [1;] propose particular Sturm-Liouville 
problem (on a finite interval) which is obtained using a 
specific coordinate transformation, see Section II. A 
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As a consistency check of this method we use it to cal- 
culate the trace of the integral operator D{z, z'), defined 
in (32) of [1], and compare it with result of a straight- 
forward alternative calculation. Expressing the trace of 
D{z, z') as a sum of its eigenvalues /i; and using (33) [ij] 
we obtain 



Tr7^(z,z')-E'^'=ET 



1 



(1) 



For the limit of piecewise constant profile being ap- 
proached, ai is given in (66) [l| 

2eA(w,k||) 

a; =Q;(u;,fc||,£) = 1 -f , (2) 



Then one easily finds 



TrD(z,z') = E 



2e 



1=1 



P7r2 + (2eK) 



Z27r2 + (2eK)2 



2eK coth(2eK) — 1 



(3) 

On the other hand, it is straightforward to calculate 
the trace of D without appealing to any eigenvalue prob- 
lem 

.1/2 

TyD{z,z') = / dzD{z,z) = 

J -1/2 

= / dxD{x,x)a^{x) ^ —. (4) 

Evidently, ^ contradicts to ([3]). At the same time it 
is straightforward to prove the equivalence of the trace 
definitions used in ([3]) and ^ taking into account that 
ipa (33) [l| constitute a complete set of eigenf unctions 
of Sturm-Liouville problem (39) Moreover, it must 
be emphasized here that @ is valid independently of 
any particular profile of (j^{x) provided it satisfies the 
normalization condition (7)[1]. It is due to the fact that 
D{z, z) (or equivalently D{x, x)) is position- independent. 

We see that the spectral problem used for calculation of 
the eigenvalues of D in the case of (approximated) piece- 
wise constant profile is either ill posed, or not equivalent 
to the original problem. It gives a wrong answer ([3]) for 
the trace oi D. In virtue of (31) [l| the same must also 
hold for the trace of the operator In AC which defines the 
Casimir energy (30) [l|. Hence, it is not legitimate to use 
the eigenvalues a; ([2]) for the (correct) calculation of the 
Casimir energy (43) [l|. 

Thus, we must conclude that the calculation method 
proposed in [l| contains an internal inconsistency when 
applied to (approximations of) piecewise constant profile 
of the mirror. One can also check that an expansion of 
(67) 11] in a power series in constant A (acknowledged in 
Section III.C [l|) does not coincide with the usual per- 
turbation theory. 

These arguments unambiguously show that the result 
for the Casimir energy of a single slab (68) [Ij is presum- 
ably incorrect. 



II. ALTERNATIVE APPROACH 

In we presented a detailed treatment of the similar 
system of piecewise constant profile without reference to 
more general cases. We restricted ourselves to consider- 
ation of homogenous and isotropic infinite plane slab of 
thickness e, placed in the xiX2 plane. 

Casimir interaction of (multi) layered systems has been 
actively studied for dielectric materials (e.g., within sur- 
face modes formalism 11], or macroscopical field oper- 
ators method ^12]). However, the mathematical formu- 
lation of these problems differ from one considered here, 
and no direct comparison of the results is possible. 

For modeling of the interaction of massive scalar field 
with volume defects we followed the Symanzik approach 
0. The complete action of the model is following 



S= So + Si + J<f> 
1 
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5o = ^ / d^x{d^(bix)d^(bix)+m^^''{x)), 



Si^^J d^x9{e,X3)c^^{x) 



The distribution function d{e,X3) represents the piece- 
wise constant profile being equal to 1/e when jxa] < e/2, 
and zero otherwise. Then there is no need for introduc- 
tion of implicit variables' change defined by (21) [ij, and 
we can proceed explicitly in Cartesian coordinates. 

We consider the generating functional of Green's func- 
tions and similarly to \§\ introduce auxiliary fields defined 
on the support of the defect (this line is also followed in 
[ij] with minor generalizations). These fields satisfy free 
boundary conditions on the edges of the layer (in other 
words there are no constrains imposed on the fields). Fol- 
lowing this approach one introduces integral operators 
acting on the auxiliary fields with support in the finite 
interval, such as A^^ in (12) of [l|. In our case it reduces 
to a unity operator and there is no need for any special 
treatment of its boundary conditions. 

Functional integral takes then explicitly gaussian form 
and we are able to derive both the Casimir energy and 
the propagator, arriving at 



-1/2 



Z[J] ^ (DetQ) exp 
S = D -T^iAO)Q-^iOA), Q = l 



{ij^j}. 



(5) 



here rj ~ Ae^^. The definitions of S and Q must be 
understood in terms of integral operators. With A = 
{—d^ + m^)^^ we denote the standard free propagator of 
scalar field, and the projecting operator O acts as 



^/2 

dx I dx3tp(j>, 

-e/2 



The unity operator 1 is also defined on the defect only. 
We shall note here that our operator Q is the direct ana- 
log of the /C (24) 1]. However, the trace calculation which 
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defines the Casimir energy and is presented below, differs 
significantly in our approach. 

The Casimir energy density per unit area of the layer 
is given by 



■Trln[Q(p; 2:3,^3)], 



(6) 



where the Fourier transformation of the coordinates par- 
allel to the defect (i.e. x = {xo,xi,X2)) was performed, 

P = {P0,Pl,P2)- 

Introducing operator U = — 1 and using the defini- 
tion of Q we can express the 77-derivative of the integrand 
of ([6]) in the following form 

d„ InQ = - — . 

V 

For explicit calculation of U we note that it is pro- 
portional to a Green function of an ordinary differential 
operator: 



KpU = -77 



(7) 



& 



Kv{x,y)= [-— + V']S{x-y). 



Then one can derive the coefficients a, b as 



2p ' 



2p 



(9) 



^^{e'^P^E + pY-iE-pf 



For the energy density we have 



,3-d 



cfip n dr] 



■TtU. 



(10) 



where we introduced dimensional regularization to free 
oneself from ultraviolet divergencies (d = 3 corresponds 
to removing of regularization) , and an auxiliary mass pa- 
rameter fj,. We have chosen the lower limit of integra- 
tion over 77 to satisfy the energy normalization condition 
£\,-i^Q = 0. It can be shown that the integral is conver- 
gent at ?7 = 0. 

The calculation of trace of the operator U is straight- 
forward 



where p — y^rj + E'^, E ~ i/p^ 4- w?. 

Employing the symmetry conditions U{x, y) = 
U{y,x) = U{—x,—y) which follow from the definition 
of U we can write 

U = — r/ — h2acosh {{x + y)p) + 2bcosh {{x — y)p) . 



2p 

From ([5|) and definition of U it follows that 
U + riO/^0{l + U) = 0. 



(8) 



T.U^ r dxU{x^x) = 2he+^-^^^^^Pl^. (11) 

J-e/2 2p 

Next, putting Q into pT|) we can prove directly that 



^ 1 



-e(E+p} 



AEpS, 



(12) 



Thus, from (|10p and we obtain the following expres- 
sion for the Casimir energy 



3-d 



d'^p 



2(27r) 



d-l 



In 



4.Ep 



{e'P{E + pf-e-'P{E-pY 



r 



= /b2 + a? 



(13) 



It can be shown that p3)l is in full agreement with the 
usual perturbation theory 

It is easily to generalize (IT^ for non-local translation 
invariant \ = \{x — y) as considered in fTl (see equation 
(4) there). In this case the final expression for 8 (fTS]) 
remains the same provided that A is replaced there with 
\{p) - Fourier image of \{x) defined in (26) ^. In the 
massless model (with m — 0, E = \p\) the energy £ 
obtained in such a waycan be compared with p^~'^£q for 
£0 presented in (67) [i| for plane slab. One can easily 
check that these results do not coincide. 

To compare ([T^ with results of Q and Q one needs to 
replace the zeta-function regularization used there with 



dimensional one applied in our approach. Using (15) of 
[3] one can obtain a generalization of (8) [3i] (written for 
d = 3) to d-dimensional case in the form 



3-d 



2(27r)° 



Jm ok 



(14) 



where = 2tt''/'^ /T{d/2) is the volume of the {d - 
l)-dimensional sphere in the d-dimensional space, and 
sii(zfc) for considered case of plane slab is defined by 
(22) Q. Integrating by part in p4)) and changing inte- 
gration variable k — ^y\p[^ + rr^P, one can check that the 
right hand sides of and (HH) coincide (up to change 
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of notation L = e, Vq = ij). In a similar way one can 
verify that result obtained in Q for the plane slab for 
d = 2 with help of zeta-function regularization agrees 

with in]). 

We hope that the arguments given above are sufficient 
to conclude that the calculation methods proposed in [l| 
need to be thoroughly verified in order to expose and 
eliminate theirs defects. Then it may be possible to em- 
ploy effectively the basic ideas of jl| in calculation of 
Casimir energy in the models of quantum field theory 



with nontrivial background. 
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